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Abstract 



o 



Context. The luminosity-metallicity relation is one of the fundamental constraints in the study of galaxy evolution; yet none of the 
. relations available today has been universally accepted by the community. 

QLh' Aims. The present work is a first step to collect homogeneous abundances and near-infrared (NIR) luminosities for a sample of dwarf 

irregular (dIrr) galaxies, located in nearby groups. The use of NIR luminosities is intended to provide a better proxy to mass than the 
blue luminosities commonly used in the literature; in addition, selecting group members reduces the impact of uncertain distances. 
Accurate abundances are derived to assess the galaxy metallicity. 
^ ■ Methods. Optical spectra are collected for H ii regions in the dirrs, allowing the determination of oxygen abundances by means of the 

' temperature- sensitive method. For each dIrr galaxy //-band imaging is performed and the total magnitudes are measured via surface 

photometry. 

Results. This high-quality database allows us to build a well-defined luminosity-metallicity relation in the range -13 > Mh > 
^ ■ -20. The scatter around its linear fit is reduced to 0.11 dex, the lowest of all relations currently available. There might exist a 

QQ ' diff'erence between the relation for dIrrs and the relation for giant galaxies, although a firm conclusion should await direct abundance 

determinations for a significant sample of massive galaxies. 
\^ . Conclusions. This new dataset provides an improved luminosity-metallicity relation, based on a standard NIR band, for dwarf star- 

■ forming galaxies. The relation can now be compared with some confidence to the predictions of models of galaxy evolution. Exciting 

• ' follow-ups of this work are (a) exploring groups with higher densities, (b) exploring nearby galaxy clusters to probe environmental 

I eff'ects on the luminosity-metallicity relation, and (c) deriving direct oxygen abundances in the central regions of star-forming giant 

, galaxies, to settle the question of a possible dichotomy between the chemical evolution of dwarfs and that of massive galaxies. 

OO . 
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1. Introduction in dwarf galaxies. Indeed, a well-defined L-Z relation is ob- 

H I served in spheroidal galaxies (e.g. Caldwell et al. 1 19921 ). mostly 

_Cd , According to standard scenarios of galactic chemical evolution, galaxies, in which old stellar populations dominate, and 

a luminosity-metallicity (L-Z) relation is established through ^r which the chemical evolution has stopped. Do we expect a 

galactic winds - induced by supernova explosions - removing ^^^^^^ behaviour in the case of dwarf irregular galaxies (dIrr)? 

the interstellar medium (ISM) before it has been totally con- ^he situation is more complex. First, dIrrs are still active (in 

verted into stars (seethe seminal work ofLarson 1974, andalso (he sense of star formation), exhibiting substantial gas fractions 

Tinsley & Larson[T979l Dekel & Silk[T986l Lynden-Bell[T992 and a broad range of star-formation rates (SFR). So, one would 

among others). As this process is more efficient in low-mass g^pect to encounter in Irrs, ISMs of different chemical maturi- 

(low escape velocity) galaxies, less metals should be trapped jies and rapidly evolving luminosities: all factors which would 

loosen any established relation between mass and metallicity. 

Send offprint requests to: isaviane@eso.org Second, the available sample of dIrrs is made of galaxies be- 

* Based on data collected at the European Southern Observatory, longing to groups rather than distributed in the field. It is well 

La Silla, Chile Proposal N. 70.B-0424(A,B); Based on observations ^^^^ ^^^^ jj^^ evolution of a galaxy gas content depends on its 

made at Lick (UCO) Observatory; Based on observations made with . ^j^u^ ^•••u^j^u ui 

xurn- TT u 1 1 .J • 1 J r. 1 u environment, and that gas stripping is boosted through galaxy 

the William Herschel Telescope operated on the island ot La Palma by . . . , ^ , f 11 o 1 ./ 

the Isaac Newton Group in the Spanish Observatorio del Roque de los interactions in dense groups and cluster sub-clumps. So, while 

Muchachos of the Instituto de Astrofi'sica de Canarias. one can see reasons why the terminal chemical evolution in dIrrs 



might depend on the galactic mass - similarly to the sample of 
gas-free spheroidal galaxies -, there are a number of extra pro- 
cesses blurring the situation. Hence, it is unclear whether dirr 
galaxies should exhibit a luminosity-metallicity relation. 

Under such a scenario dwarf galaxies lose substantial frac- 
tions of their mass in the course of their evolution, and should be 
major contributors to the chemical evolution of the inter-galactic 
medium (IGM). In particular, mass loss through galactic winds 
should occur in dIrr galaxies: their study contributes to probe the 
fraction of the IGM enrichment which is due to dwarf galaxies 
in general (e.g., Garnett 2002) . Losses of enriched g as in d warf 
galaxies have indeed been observed: Martin et al. (120021) de- 
tected a galactic wind of ~ 6 x 10^ At© from NGC 1569, and 
for the first time, they could measure its metallicity. The au- 
thors concluded at 3 x 10"^ At© of oxygen in the wind, almost 
as much as in the disk of the dwarf itself. Hence, the cosmic 
chemical evolution and the existence of an L-Z relation appear 
to be closely inter-related: this was an additional motivation to 
pursue the study presented here. 

The case of another dIrr, SagDIG, further illustrates this 
point. In one of our previous studies (Saviane et al. 12002b we 
measured a very low oxygen content (7.26 < 12 -\- log (0/H) < 
7.50), a result which is not compatible with a closed-box evo- 
lution. According to the model, such a low abundance would 
imply a large gas mass fraction yu = mgas/(mgas + ^stars) ~ 0.97, 
whereas the observed gas mass fraction is only ji ^ 0.86. It is 
easy to compute that 1.5 x 10^ At© of gas are missing, and, 
since this mass is smaller than that of the NGC 1569 galactic 
wind, it is plausible to conclude that SagDIG lost some of its gas 
into the IGM. 

On the other hand, a cautionary remark should be added, 
since a general consensus on the role of galactic winds in the 
evolution of dwarf galaxies and the IGM has not been reached 
yet. For example by means of chemo-photometric models Calura 
& Matteucci (120061) conclude that dIrr galaxies play a negligible 
role in the enrichment of the IGM, and the numerical models of 
Silich and Tenorio-Tagle (1998) show that that galactic winds 
never reach the escape velocity of these dwarfs. And in the case 
of dwarf spheroidal (dSph) galaxies, the 3D hydrodynamic sim- 
ulations of Marcolini et al. (120061) show that no galactic winds 
develop in these objects. 

In our attempt to place SagDIG on previously established L- 
Z relations for dIrr galaxies, we realized that the existence of 
such a relation was controversial. Some studies had concluded 
at very well-defined correlations (Skillman et al. 1989: Richer & 
McCall ll995[ Pilvugin l200T]) . others have found only mi ld rela- 
tions with substantial scatter (e.g. Skillman et al. I2003bl) . or no 
correlation at al (Hidalgo Gamez & Olofsson 119981 Hunter & 
Hoff^man ll999l) . Perhaps the major source of confusion in these 
investigations comes from the ill-defined samples used in the 
analysis. Abundance data are taken from diff'erent sources (with 
spectra of variable quality and processed through diff'erent re- 
duction and analysis techniques); apparent luminosities and dis- 
tances are taken from catalogs, meaning that they are largely 
approximate. Moreover, although the L-Z relation is expected 
to depend on the environment, this is rarely taken into account 
in the analysis. Another main source of uncertainty is the com- 
mon use of the blue absolute magnitude as a tracer of the mass. 
Optical luminosities can in fact be extremely misleading in Irr 
galaxies, because of their star-bursting activity (e.g. Tosi et al. 
fT992l or Schmidt et al. lT995]) . Already, Bruzual & Chariot (IT983]) 
showed that at 400 nm a 1 Myr old burst is ~ three orders of mag- 
nitude brighter than an underlying old (15 Gyr) stellar popula- 
tion of comparable mass. In other words, a dwarf galaxy host- 



ing a recent starburst could be as luminous, in the blue, as a 
galaxy orders of magnitudes more massive but lacking a recent 
star formation episode. In comparison, a recent starburst is at 
most 10 times brighter in the near-infrared (NIR) than its under- 
lying old population of same mass. Therefore, the NIR window 
being more stable with regard to star formation episodes is more 
appropriate for probing the galaxy basic properties, such as its 
mass. 

In order to address these issues, we embarked on a medium- 
term project aimed at gathering nebular oxygen abundances and 
NIR luminosities for a sample of dirrs belonging to the three 
nearest groups of galaxies. In this way we can test the exis- 
tence of an L-Z relation in well-defined environments (charac- 
terized by the group density), for which the scatter in apparent 
distance modulus is low, and for which an homogeneous set of 
abundances can be obtained. 

Since we started the project, a few investigations have been 
published, in relation with the context presented above. The 
first L-Z relation using NIR luminosities (aside the work by 
Saviane et al. l2004l) is that by Salzer et al. (120051 hereafter S05). 
Their sample is dominated by massive galaxies, and their rela- 
tion was derived using 2MASS data for the luminosities, and 
proprietary spectra from the KPNO International Spectroscopic 
Survey (KISS) for the metallicities. The 2MASS survey having a 
relatively shallow limiting magnitude, for a few additional dwarf 
galaxies NIR photometry was supplemented by the authors. At 
the other extreme, the sample assembled by Mendes de Oliveira 
et al. (120061) is entirely made of dwarf galaxies. They gathered K- 
band luminosities and metallicities - obtained through the direct 
method - for 29 dIrrs. T he NIR data are from 2MASS or from 
Vaduvescu et al. (120051) . and metallicities have b een picked up 
from a variety of sources. Finally Lee et al. (12006b have used the 
Spitzer Infrared Array Camera to compute 4.5 yum luminosities 
for ~ 30 nearby dIrrs, the distance of which have been derived 
using standard candles. With oxygen abundances collected from 
the literature, they constructed a 4.5 jim L-Z relation. 

Although these studies have certainly improved the situation, 
they do not represent the ideal case yet. S05 abundances are de- 
rived through a new but indirect method, and the sample is very 
scarce in dwarf galaxies. Yet, S05 is the only study which in- 
cludes //-band photometry, so a comparison with their L-Z rela- 
tion is carried out later in this paper, in Sect. 15.31 As extensively 
discussed in Sect. 15.21 the use of heterogeneous data from the 
literature may produce a diff'erent L-Z relation than the one ob- 
tained with an homogeneous dataset, a likely consequence of 
distance uncertainties. Therefore, we anticipate that the relation 
by Mendes de Oliveira et al. (12006b will need to be revised once 
a better controlled sample, with ^-band photometry, becomes 
available. Finally, the use of a non-standard passband d oes not 
allow an easy comparison of the L-Z relation by Lee et al. (12006b 
to other L-Z relations. For example they need to make several 
successive assumptions in order to convert their luminosities into 
masses, and then compare the mass-metallicity relation to the 
SDSS one. For all these reasons, we believe that so far our ap- 
proach stands as the one with the smallest number of uncertain- 
ties and the broadest applicability. 

The paper is structured as described below. In Sect. 2 we ex- 
plain how the targets were selected, and give a brief account of 
the data reduction techniques - described in more detail in the 
appendices -. The computation of chemical abundances is de- 
scribed in Sect. 3, where we also compare our results with pre- 
vious abundance determinations. The resulting NIR luminosity- 
metallicity relation for dIrr galaxies is presented in Sect. 4. It is 
discussed in Sect. 5, with a possible explanation of its origin pre- 



Table 1. Target list 



Main ID 


Other IDs 




M81 group 




DDO 42 


NGC 2366 


UGC 03851 


KDG 52 


M81 dw A 


PGC 023521 


DDO 53 


UGC 4459 


VII Zw 238 


UGC 4483 






KDG 54 


UGC 4945 


UGCA 158 


candidate 5 


[FM2000] 5 




BK 3N 


PGC 028529 




DDO 66 


UGC 05336 


Holmberg IX 


DDO 82 


UGC 5692 




DDO 165 


UGC 8201 


VII Zw 499 




Sculptor group 




ESO 347-G017 






UGCA 442 


ESO 471- G 006 


AM 2341-321 


ESO 348-G009 


AM 2346-380 


UKS 2346-380 


NGC 59 


ESO 539- G 004 




ESO 294-GOlO 


AM 0024-420 




ESO 473-G024 


AM 0028-23 




AM0106-382 


LEDA 166061 


[KK98] Oil 


NGC 625 


ESO 297- G 005 


AM 0132-414 


ESO 245-G005 


AM 0142-435 


A143 



sented in Sect. 5.1 and its comparison with a relation based on 
literature data in Sect. 5.2. In Sect. 5.3 we discuss the possible 
dwarf vs. giant galaxy dichotom y, thro ugh a comparison of our 
results to those by Salzer et al. (12005b . Finally a summary and 
the conclusions of this study are provided in Sect. 6. 

2. The data set 

2.1. Target selection 

The next two sections describe how the targets were selected 
in each of the two groups. The most popular galaxy names are 
listed in Table [T] together with other aliases. 

2.1.1. Sculptor group 

Sculp tor gro up Irr galaxies were selected from the list of Cote 
et al. J1997t hereafter C97). The seven galaxies with Ha detec- 
tion and radial velocities Vr < 1000 km s"\ were included in 
our target list shown in Table [21 Note that the C97 catalog con- 
tains another set of galaxies with velocities between 1000 and 
2000 km s"^ as well as sources with velocities up to ~ 17000 km 
s"^ In addition to our main list, we included ESO 294-GOlO, 
which has been classified as intermediate dSO/Im by Jerjen et 
al. (119981 hereafter J98) and, although C97 measured an Ha ve- 
locity Vr = 4450 km s"^ a much smaller velocity is quoted in 
NED. Moreover, J98 have detected [O iii] and Ha emission, al- 
beit very weak. Yet, as no emission was detected during our Ha 
pre-imaging, we did not carry over further observations: it is pos- 
sible that J98 detected a PN with Ha emission below our detec- 
tion threshold. Bouchard et al. (12005b re-classified this galaxy as 
dSph/dlrr. 

The main observational campaign was carried out in October 
2002: NIR imaging of the seven targets was collected, while bad 
weather undermined the spectroscopic observing run and spec- 
tra could be collected for five targets only. A few months af- 
ter the completion of our observations, Skillman et al. (12003 at 
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Figure 1. Erom left to right, we show the morphology of the 
Sculptor group targets in the optical (DSS), infrared H (SOEI), 
and a false-color montage of Ha and R (EE0SC2) images. The 
reddish, diff'use nebulae in the right image are H ii regions, and 
the vertical ticks mark those for which spectra have been ob- 
tained. The galaxy luminosity decreases from top to bottom pan- 
els, with the exception of ESO 294 GO 10, which was not ob- 
served in the NIR. In all images. North is up and East to the left, 
and the field of view is ~ 5' x 5'. 



hereafter S03a) identified other dwarfs showing Ha emission 
despite a null detection by C97, bringing the total number of 
actively star-forming dirrs in the Sc ulptor gr oup up to twelve. In 
a companion paper, Skillman et al. (I2003b[ hereafter S03b) de- 
rived oxygen abundances for five of these galaxies (reported in 
Table [7]), two of them being in common with our original sam- 
ple. So, to improve the comparison with S03b, in August 2003 
we complemented our data set with NIR imaging and long- slit 
spectroscopy for ESO 473-G024. Einally, ESO 245-G005 was 
re-observed in October 2003, in order to measure an H ii region 
closer to the center of the galaxy. 



Table 2. Observations of Sculptor group targets. Additional data are the radial velocity from Cote et al. (11997b , and the radial 
velocity from the NASA extragalactic database. No emission was detected in the Ha pre-imaging for ESO 294-GOlO. 





Target list 








Journal of the observations 






Galaxy 


RA. 


Dec. 


Vr [km/sec] 


EF0SC2 


SOFI 






J2000 


C97 


NED 


Date 


t [s] 


Date 


t[s] 


airm. 


ESO 347-G017 


23:26:56.0 


-37:20:49 


659 


692+4 


13 Oct 2002 


3 X 1800 


15 Oct 2002 


3420 


1.10 












14 Oct 2002 


2x 1800 


16 Oct 2002 


300 


1.02 


UGCA 442 


23:43:45.5 


-31:57:24 


283 


267+2 






16 Oct 2002 


3600 


1.08 


ESO 348-G009 


23:49:23.5 


-37:46:19 


628 


648+4 






16 Oct 2002 


3600 


1.09 


NGC 59 


00:15:25.1 


-21:26:40 


357 


382+60 


13 Oct 2002 


4x 1800 


15 Oct 2002 


4080 


1.05 
















16 Oct 2002 


300 


1.02 


ESO294-G010 


00:26:33.4 


-41:51:19 


4450 


117+5 












ESO 473-G024 


00:31:23.1 


-22:46:02 




541 + 1 


19 Aug 2003 


3 X 1200 


10 Aug 2003 


1560 


1.10 


AM 0106-382 


01:08:21.9 


-38:12:34 


645 


645+10 


14 Oct 2002 


4 X 1800 


15 Oct 2002 


3480 


1.26 
















16 Oct 2002 


1800 


1.15 


NGC 625 


01:35:04.2 


-41:26:15 


415 


396+1 


14 Oct 2002 


3x900 


15 Oct 2002 


1020 


1.71 
















16 Oct 2002 


900 


1.78 


ESO 245-G005 


01:45:03.7 


-43:35:53 


389 


391+2 


13 Oct 2002 


2 X 1800 


16 Oct 2002 


2460 


1.32 



14 Oct 2002 3 x 1800 
14 Oct 2003 3 X 1200 



DSS H Ua^ 



HI 


* 






♦ ■ » 

♦ 




■'t 









Figure 1. (continued) 



In summary, with the observations presented here there are 
now eight Sculptor group dwarfs with measured abundances, of 
which three are in common between our study and S03b. NIR 
data are available for all eight objects. 



2.1.2. M81 group 

A first list of M81 group members has been compiled by 
Karachentseva et al. (119851 hereafter K85), by inspecting photo- 
graphic plates taken at the 6-m telescope. It includes 14 probable 
and 1 1 questionable members (on the basis of radial velocity and 
visual appearanc e), all o f them of dirr type. Later, Karachentseva 
& Karachentsev (119981) added six new possible members in the 
dE class, by searching POSS-II and ESO/SERC films. A search 
based on CCD imaging was carried out by Caldwell et al. (11998b , 
who, however, did not publish their catalog. Finally, Eroebrich & 
Meusinger (120001 hereafter EMOO) discovered six more candi- 
dates in a survey that employed digitally stacked Schmidt plates. 
The candidates are equally divided into dE and dIrr types. A new 
version of the catalog of M8 1 member candidates was recently 
published by Karachentsev et al. (120021 hereafter K02), show- 
ing a few additions and a few deletions when compared to the 
previous catalog. 

Ten M81 group dwarfs have been observed from 2001 to 
2002 (Table [3]). Our first run was with the Kast spectrograph 
at Lick: the instrument off'ers limited imaging capabilities, so we 
selected the target H ii regions with the following strategy. The 
slit was first placed on the brightest nebulosity within a target 
dwarf and a short-exposure spectrum was taken. If it showed a 
conspicuous Ha line, then the integration was continued, other- 
wise we moved to another candidate region and tried again. If no 
emission was detected in two to three regions (depending on the 
galaxy luminosity), we moved to the next dwarf target. With this 
procedure, emission was detected in only half of the galaxies, 
namely DDO 42, DDO 53, UGC 4483, DDO 82, and KDG 54 
(which turned out to be a background, higher redshift object). No 
emission was detected in KDG 52, BK 3N, DDO 66, DDO 165, 
nor in the candidate 5 of Eroebrich & Meusinger (12000b . For 
objects which did show emission lines, the spectroscopic obser- 
vations were followed-up by NIR imaging performed with the 
INGRID camera at the ING/WHT telescope on La Palma. 

After the completion of the spectroscopic observations, three 
Ha surveys including M81 group galaxies were published by 
Hunt er & El megreen (l2004l) . Gil de Paz et al. ((20031) . and James 
et al. (120041) . These surveys have detected emission-line regions 
in, among others, DDO 42, DDO 53, UGC 4483, DDO 66, 
DDO 82, and DDO 165. So they confirm our detections, but they 



Table 3. Observations of M81 group targets. The group membership has been estabUshed in K85, except for the candidate 5 
which comes from Froebrich & Meusinger (2000). Additional data are the radial velocity from the NASA extragalactic database. 
Exposure times for the Kast spectrograph observations are flagged with a code indicating the instrument setup: B=blue arm with 
grism 830/3460; R=red arm with grating 300/4230; BR=both arms with dichroic D46; F=red arm without dichroic; HR=red arm 
with grating 1200/5000 centered around Ha. A 'n.e.' in the INGRID observation means that no emission was detected and so, no 
NIR imaging had been performed. 





Target list 






Journal of the observations 




Galaxy 


R.A. 


Dec. 


Vr 




Kast 


IJN LrKiJJ 




Main TD 


J2000 


ricm/<sPpl 




t[s] 


Date t [s] 


airm 


DDO 42 


07:28:53.7 


+69:12:54 


+80 


17 Dec 2001 

1 7 Der 9001 

1 7 Hpp 9001 

14 Tan 9009 
IH- Jail z.\J\JjL 

14 Tan 9009 

1 ^ Tan 9009 
1 J Jail jL\J\JjL 


600+2x1200 RB 
1800 B + 1800 HR 
600 F + 600 HR 
600+2x1200 RB 
1800 B + 1800 HR 
1800 B + 1800 HR 


1j Jan zUUJ JoUU 
lo Jan zUUJ JUU 


1.36 
1.46 


KDG 52 


08:23:56.0 


+71:01:45 


+ 113 


17 Dec 2001 


2 x300 RB 


n.e. 




LJUkJ od 


08:34:07.2 


+66:10:54 


-1-90 


1 7 T^pp 9001 
17 Dec 2001 


J A i OUU IN. O 

1800 HR+ 1800 B 


1 f\ Tan 900^ l^f\C\ 
iU Jail zuvjj) ZiUVJ 


1 


UGC 4483 


08:37:03.0 


+69:46:50 


+ 156 


14 Jan 2002 
14 Jan 2002 


2X1200RB 

1800 B + 1800 HR 


15 Jan 2003 960 

16 Jan 2003 3210 


1.29 
1.34 


KDG 54 


09:22:25.2 


+75:45:57 


+659 


14 Jan 2002 


2x(1200+1800)RB 


higher redshift 




candidate 5 


09:39:02.0 


+69:25:01 








n.e. 




BK3N 


09:53:48.5 


+68:58:08 


-40 






n.e. 




DDO 66 


09:57:30.1 


+69:02:52 


+46 






n.e. 




DDO 82 


10:30:34.8 


+70:37:14 


+56 


15 Jan 2002 
15 Jan 2002 


300 + 2x600 RB 
3x1200 RB 


16 Jan 2003 3600 


1.41 


DDO 165 


13:06:24.8 


+67:42:25 


+31 






n.e. 





also unveil emission-line regions in DDO 66 and DDO 165, for 
which our trial-and-error technique gave a null result. 

2.2. Optical spectrophotometry 
2.2.1 . Sculptor group 

The observations summarized in Table [2l were carried out us- 
ing EF0SC2 ( ESO F aint Object Spectrograph and Camera, 2; 
Buzzoni et al. 11984b at La Silla. Most spectra were taken on 
October 13 and 14, 2002, an additional spectrum of ES0245- 
G005 was collected on October 14, 2003, and finally ES0473- 
G024 was observed in August 2003. EF0SC2 is a multi-mode 
instrument working at the Cassegrain focus of the ESO/3.6m 
telescope, allowing both imaging and long-slit, low-resolution 
spectroscopy. A camera images the aperture onto a 2048 x 
2048 px Loral/Lesser, thinned, and UV-flooded CCD. Its pixel 
size is 15yum, which corresponds to Qf!l51 on the sky, for a total 
field of view of 5.4 x 5.4 Under normal operation, the CCD 
is read out by the left amplifier at a gain of 1.3e"/ADU, and 
since numerical saturation occurs at 2^^ADU, it remains well 
below the full well capacity of 104, 000e"px"^. The dark current 
is 7e"px"^hr"\ comparable to the readout noise of ~ 10 e"px"^ 
The CCD was binned at 2 x 2 both in imaging and spectroscopy 
mode. Indeed, rebinning allows the line profile of the calibra- 
tion arcs to be sampled with more than four pixels, even with 
the highest dispersion grisms and a Of! 5 wide slit, so for spectro- 
scopic observations there is no real advantage in using the 1 x 1 
binning mode. 

For each galaxy, the target H ii regions were identified by 
subtracting a 5 min R image from a 15 min Ha image, and the 
slit was then centered on the brightest region. Maps of the H ii 
regions are shown in Fig. [T] During the night, two spectropho- 
tometric standard s (Feige 110 and LTT 3218) from Oke ([T990b 
and Hamuy et al. (119921 [T994b were observed. Other calibration 



frames were taken in the afternoon preceding each observing 
night (bias, darks, HeAr arcs, and dome flatfields). Grism #11 
was used without an order sorting filter, yielding a dispersion of 
4.2 A px"^ and giving a range from 3400A to 7500A. The second 
order spectrum (appearing at wavelengths longer than ~ 6800A) 
is only a few percent of the first order, so the contamination is 
negligible for the target spectra. It is however important for the 
determination of the response function. Indeed, it can be seen in 
the spectrum of the extremely blue standard Feige 110 (spectral 
type DO), but it disappears already in the spectrum of LTT 3218, 
although its spectral type is DA. The latter standard was then 
used to compute the response function. Spectra of the target re- 
gions were obtained with a width slit, while the 5'' width slit, 
aligned along the parallactic angle, was used for observing the 
spectrophotometric standards, in order to minimize slit losses. 
The resolution with the V width slit is ^ 13A FWHM (from the 
HeAr lines near 5000A). All slits cover the entire EF0SC2 field 
of view. 

The reduction and analysis of the spectroscopic data was per- 
formed in a standard way, with the main steps: (a) 2D wave- 
length and distortion calibration; (b) extraction of spectra using 
the wings of Ha to define a good window; (c) flux calibration; 
(d) measurement of line fluxes in the midas/alic^ framework; (e) 
correction of fluxes for underlying absorption; (f) correction for 
internal reddening of the region using the Ha/H/? ratio; (g) com- 
putation of Te using [Out] lines /14363 and /1/14959, 5007; and 
finally (h) computation of nebular abundances using iraf/ionicQ. 



^ Eso-MiDAs is the acronym for the European Southern Observatory 
Munich Image Data Analysis System which is developed and main- 
tained by the European Southern Observatory 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Figure 2. The NIR, //-band, morphology of the M8 1 group dirrs 
observed at WHT with ingrid (right column) is compared to their 
optical appearance as seen in DSS images (left column). The 
galaxy luminosity decreases from top to bottom panels. In all 
images, North is up and East to the left, and the field of view is 
~ 6' X 6' in the left column and ~ 5' x 5' in the right column. 



In some cases the Ha line profile showed more than one re- 
gion: in this case the two brightest regions were extracted (see 
Table 131). The spectra of the brightest H ii regions that enter the 
L-Z relation for each galaxy are shown in Fig.O 

Full account of the data reduction is reported in AppendixlAl 
The reddening corrected fluxes for each useful region are listed 
in TablelH for the regions where [Oiii]/14363 has been detected, 
we list their physical parameters and abundances in Table [6j 

2.2.2. M81 group 

The spectra of the M81 group dwarfs were obtained with the 
Kast Double Spectrograph at the Cassegrain focus of the Shane 
3m Telescope, at Lick observatory. The instrument consists of 
two separate spectrographs - one optimized for the red, and the 
other for the blue -. Dichroic beam-splitters and separate CCD 
detectors allow simultaneous observation. Both arms use a UV- 
flooded Reticon 1200x400 px device, with 2.1 jim pixels cover- 
ing O'.'S on the sky; the ADC sets a maximum ADU count of 
32000, with a gain of 3.8 e"/ADU, and the read-out noise is 
6 e". The CCDs have a 40% quantum efficiency at 3200A.We 
used a V!5 wide slit. The observations were carried out simulta- 
neously with the blue and the red arms, using dichroic D46. In 




4000 5000 6000 7000 

A[A] 

Figure 3. For each dirr in the Sculptor group sample we show 
the spectrum of one of the H ii regions that enter the L-Z relation. 
The spectra are masked near the strong night sky line [O i]/15579, 
which leads to poor sky subtraction. The reddest, noisiest, part 
of the spectrum is omitted for two objects. For display purposes, 
each flux has been multiplied by 10^^ and a constant has been 
added. 
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Figure 4. For three galaxies in the MS 1 group sample the spec- 
trum of the target H II region is shown, with a label identifying 
the host galaxy. The spectra of other galaxies are of poor qual- 
ity and have not been used for abundance determinations. For 
display purposes, each flux has been multiplied by 10^^ and a 
constant has been added. 



Table 5. Reddening-corrected fluxes for aU H ii regions in M81 group galaxies. 
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Notes: (^) Sii /15640? (^) The reddening constant is computed by the expression: iHa/Injs = iHao/lHjso^O ^^f^^"^ fW)]^ where /q and / denote the 
flux before and after extinction, and f{A) is the extinction curve. 



one occasion we also obtained a spectrum with the red arm alone 
(without dichroic) to check the flux normalization in the overlap 
region of the red and blue spectra. Table[3]summarizes the obser- 
vations, carried out in December 2001 and January 2002. Due to 
poor weather conditions, and the overheads of the trial-and-error 
technique, only six galaxy targets were observed in three nights, 
of which one is the higher redshift object KDG 54. 

The reduction and analysis of the spectroscopic data is sim- 
ilar to that of EF0SC2 (see Appendix O, although more steps 
are needed to correct for the large flexures and to ensure a proper 
normalization of the spectra in the two arms (see Appendix |B]). 
Spectra with adequate S IN ratio could be obtained for three 
galaxies, they are shown in Fig.lH 

The reddening corrected fluxes for each useful region are 
listed in Table [51 and for the regions where [Oiii]/14363 was 
detected, their derived physical parameters and abundances are 
provided in Table [6l 

2.3. NIR imaging 

The imaging part of the project was carried out with the 
4m William Herschel Telescope (WHT) and the 3.5m New 
Technology Telescope (NTT). The instrument INGRID (Isaac 
Newton Group Red Imaging Device) was used at the former, 
and SOFI (IMorwood et al. 119981) at the latter. Both instruments 
make use of 1024x1024 HgCdTe Hawaii arrays, and the chosen 
scales were respectively of 0^.^238 px"^ and 0^.^288 px"^. 



The observations were performed with the typical pattern for 
infrared imaging: iterating every 1-2 minutes between the ob- 
ject and a nearby patch of clear sky, usually 3-5' away from the 
galaxy. Each image was averaged over a number of short integra- 
tions, usually 10-30 sec, to avoid saturation of the array from the 
combination of the sky background and any bright foreground 
objects. Small ofl'sets - of order of 5-10''- were introduced be- 
tween successive object or sky images to ensure that the target 
was placed on difl'erent locations of the array. This helps remov- 
ing cosmetic defects and improves the flat-fielding. A log of the 
observations is presented in Tableland in Table [5] 



Detailed account of the data reduction and calibration is pro- 
vided in Appendix O The main steps of the NIR surface pho- 
tometry are: (a) accurate flat-fielding with polynomial fitting of 
the background; (b) masking of foreground stars; (c) deriving the 
light-profile by integration in elliptical annuli. Errors on the to- 
tal magnitudes have been estimated by changing the background 
level by ±lcr, integrating the two new growth curves, and then 
computing the diff'erence with the magnitude obtained from the 
average growth curve. In turn, the r.m.s. fluctuation of the back- 
ground was estimated by measuring the background itself in 10- 
20 independent areas. The total magnitudes of all galaxies are 
listed in Table [71 except for DDO 82, dropped out because no 
estimate of its abundance could be obtained; its magnitude is 
H= 10.83 ± 0.37 mag. 



Table 6. Physical parameters and abundances for Hii regions in the dirrs of the Sculptor and M81 groups, for which the electron 
temperature could be determined. The and S^/H^ abundances are the average of the two values obtained from the two lines 

of the doublet separately (AA5001, 4959, and AA6116, 6731, respectively) 
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2.84 + 0.18 


3.82 + 0.53 


(xlO^) 


1.90 + 0.96 


8.67 + 7.45 


3.54 + 1.03 


1.06 + 0.47 


0.31+0.06 


0/H (xlO^) 


4.17 + 0.96 
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log(S/0) 


-1.89 + 0.23 
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n{S II) (cm-3) 
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135 + 63 
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8.37 + 0.07 
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8.39 + 0.08 


8.42 + 0.12 
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-1.29 + 0.41 
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153 + 98 




156 + 139 






ICF 


1.43 + 0.37 




1.95 + 0.86 






S/H (xlO^) 


23.06 + 9.87 




31.80 + 14.02 






12 + log (S/H) 


7.32 + 0.20 




7.46 + 0.21 






S/O(xl02) 


0.97 + 0.58 




1.27 + 0.78 






log (S/0) 


-2.01 + 0.20 




-1.90 + 0.21 







3. Derivation of the chemical abundances 

Since the main purpose of this study is to use oxygen abundances 
in the search for a correlation with the galaxy luminosity, we 
defer a complete discussion of the elements other than oxygen 
to a future paper. Yet, the abundances of two other important 
elements (nitrogen and sulfur) have been derived as described 
in this section, and are listed in Table [S] The abundances have 
been computed with the so-called direct method (e.g. Osterbrock 
fT989l) . 

The O III region electron temperatures Tq, electron densities 
Hq, and abundances have been computed using tasks within the 
NEBULAR package of iraf. Since temperature and density cannot 
be derived independently, we start with suitable initial values 
and iterate the task temden until the values stabilize within the 
errors: the temperature is computed using the [Oiii] line ratio 
(^4959+.15007)/.14363, and the density is derived from the [S ii] 
line ratio /16716//16731. Using these values for Tq and n^, for 
every ion the abundance is computed with the ionic task, with 
central wavelengths, line ratios, and errors listed in Tableland 
Table [5] for the Sculptor and M81 groups, respectively. The tol- 
erance on the central wavelength is generally taken as lA, ex- 
cept in cases where the line is actually an unresolved doublet 



or blend. This happens for the following multiplets: [Oi]/16364 
(tolerance=2A), [Oii]^3727 (4A), [Oii]^7325 blend (lOA), and 
the [S ii]^4072 blend (7 A). 

The total oxygen abundance is computed as 
0/H = 0++/H+ + 0+/H+, i.e. neglecting the usually small 
contribution from O^^ (Skillman & Kennicutt 119931) . On the 
other hand, higher ionization stages of sulfur do contribute 
to the total abundance (Dufour et al. 11988b , so an ionization 
correction factor (ICF) has to be computed and taken into 
account. This is extensively discussed in Garnett (11989b , who 
computed a series of photoionization models and compared 
them with the relati on (S+ + S++)/S = [1 - (0+/0)^]i/^^ro- 
posed by Stasinska (11978b with a = 3 and French (1981) 
with a = 2. Neither choice of a gives an entirely satisfactory 
result, however we have found that for a = 2.6 the relation 
reproduces Garnett's models in the range 0^/0 > 0.2, which 
corresponds to the case of all our Hii regions. Thus, sulfur 
abundances are computed as S/H = ICF x (S^ + S^^), with 
ICF = [1 - (1 - 0+/0)2-6]-i/2-6; the correction factors are Hsted 
in Table [6l Finally nitrogen abundances are computed assuming 
(N/0) = (N^/0^), and then using the nitrogen to oxygen ratio 
in the equation (N/H) = (N/0) x (0/H). 
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Figure 5. Comparison of our abundance determinations to liter- 
ature values, for H ii regions in galaxies with independent - and 
direct - determinations of 12 + log (0/H). The dotted line marks 
the median deviation of -0.1 dex. The galaxy name is followed 
by the reference acronym, and its subscript is the region number 
in the reference paper. 



3.1. Comparison with previous abundance determinations 

A summary of our measurements is given in Table Ul and, for 
each galaxy, the best literature direct measurement of the oxygen 
abundance is reported as well. If no direct measurement exists, 
then we list the best indirect measurement (i.e. the temperature 
cannot be estimated, so abundance estimates rely upon semi- 
empirical relations making use of different line ratios). There 
are five galaxies in our sample for which the abundances can be 
compared to values already published: this is shown in Fig. [5] 
There is a general agreement except for DDO 42, which exhibits 
an internal metallicity gradient and is discussed in more detail 
in Sect. [H Excluding this galaxy, the dispersion in the value dif- 
ferences is only 0.06 dex. However, the median of the litera- 
ture values is 0.1 dex lower than ours (dotted line in the figure). 
We tried to understand this discrepancy by examining NGC 625, 
which shows the largest diff'erence. As a first check, we com- 
pared our fluxes to those of S03b for the same nebula in this 
galaxy: the two studies are fairly consistent, our fluxes being 
only 3% smaller on average and with a dispersion of 14%, and 
no trend with wavelength. The exception is the [S iii] /16312 line, 
for which our flux is ~ 2 times larger than that quoted in S03b. 
With respect to the oxygen lines, our fluxes are 17% and 25% 
smaller for the [On] A3121 and [Oiii] /14363 lines, respectively, 
while the [Oiii] /14959 and /15007 line fluxes are only 7% and 
6% smaller, respectively. Such flux discrepancies are entirely 
within the errors we quote, yet they are responsible for most 
of the abundance difl'erences. Indeed, using the fluxes listed in 
S03b we obtain 12 -h log(0/H)= 8.26 dex for NGC 625 (region 
5), so the difl'erence is reduced from 0.18 dex to 0.07 dex, i.e. it 
is now within our Icr error on the abundance. 

Since the measured green [O iii] line fluxes are very close to 
those of S03b, we discuss only the blue [Oiii] /14363 and [On] 
lines here. The difl'erence in flux for the [O ii] A3121 line could 
be due to the fact that for A < 4000A the EF0SC2 response 
function starts dropping steadily, and large corrections have to 
be applied. On the other hand, comparing the response functions 
obtained at difl'erent epochs shows that the uncertainty on our 
response curve is at the level of a few 0.1%. S03b do not show 
their response function, but since they are using an order sorting 
filter with 50% cut at 3600A it is very likely that their correction 
to the [O ii] line is quite large as well (almost a factor of two). In 
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Figure 6. Metallicity of Hii regions vs. absolute H luminosity 
of galaxies in the Sculptor and M81 groups. The upper panel 
shows all direct oxygen measurements listed in Table [71 while 
in the lower panel only selected regions are plotted with filled 
diamonds (see text). The dashed line is a weighted least-squares 
fit to the data, and the dotted lines feature the ±lcr dispersion 
around the fit. The galaxy IDs in the upper panel are followed 
by the acronym of the paper where the abundances have been 
published (see Table [7]), and the subscript is the region number 
in that paper. To separate the labels, small random vertical off'sets 

have been introduced. The open circles at Mh 15.5 represent 

the two components of I Zw 1 8 for the two distances listed in 
Table |7l while the open square represents SBS 0335-052. 



any case, adopting the S03b [O ii] line flux leads only to a minor 
increase in abundance of 0.03 dex. Most of the abundance de- 
pendence is then in the [O iii] /14363 line, and indeed changing its 
flux and adopting the S03b value we get 12 -h log(0/H) = 8.19. 
In this case, the flux quoted in S03b might be more reliable than 
ours. Indeed, we have to deblend this line from Hy, while on the 
spectra published by S03b the line is well separated from Hy. 
Since [Oiii] /14363 becomes weaker as the oxygen abundance 
increases - hence the temperature decreases - this could explain 
why the most deviant abundance value is that corresponding to 
NGC 625, the most metal-rich galaxy in our sample. 

4. The NIR luminosity-metallicity relation 

Table[7]summarizes the galaxy parameters which have been used 
to create the L-Z relation. Besides oxygen abundances and NIR 
luminosities, the table reports on distances an d total extinctions 
in the H band, taken from Schlegel et al. (11998b . Distances 
based on standard candles are available for only a fraction of 
the targets, so in most cases we have used approximate val- 
ues. In particular for the Sculptor group, distances were taken 
from S03a (computed with the method of J98). This method 
relies on the recession velocity, and is claimed to provide dis- 
tances within 20% of the direct measurement. The three galaxies 
with distances computed with the red-giant branch tip (RGBT) 



Table 7. Data for the L-Z relation of Sculptor and M81 dirrs, taken both from the literature and from the present work. The 
abundances, computed with the method identified in the third column, have been extracted from the publication identified by its 
acronym in the fourth column (see the bibliography); SOS is this paper. The abundance methods are D=direct, M=McGaugh, 
P=Pilyugin, and I=indirect. Notes on abundances are: (a) average of 3 regions; (b) average of two methods; (c) central region; (d) 
external region. Distances have been computed with the method identified in the sixth column, and have been extracted from the 
publication identified by its acronym in the seventh column (see the bibliography). The distance methods are J98, assuming an error 
of 20% and RGBT - red g iant branch tip -. Total apparent //-band magnitudes are from this work, and the //-band extinctions are 
from Schlegel et al. ([T998b . 



Name 


iz-l-10g(U/rlj 


Method 


Kei 


JJ [IVlpCj 


Method 


Kei 


Ah 


WH,tot 


Mh 








Sculptor 


group 












ESO 347-G017 


o f\f\ . r\ r\c 

8.00 ± 0.05 


D 


S08#l 


6.99 ± 1.40 


J98 


S03a 


0.010 


11.79 ± 0.50 


-17.40 ± 0.66 




/.Q/ ± 0.20 


P 


S08#l 
















7.90 ± 0.09 


D 


S03b#5 














UCjCA 442 






S08 


4.27 + 0.52 


RGBT 


K03 


0.010 


1 1 '~if\ 1 rv Cf\ 

11.70 ± 0.50 


-16.45 ± 0.56 




7.72 ± 0.03 


D 


S03b 














ESO 348-G009 






S08 


6.52 ± 1.30 


J98 


S03a 


0.008 


12.80 ± 0.61 


-16.24 ± 0.75 




7.89 ± 0.03 


M 


S03b, L03 
















8.11 ± 0.05 


P 


S03b, L03 














NGC 59 


8.39 ± 0.08 


D 


S08#l 


4.39 ± 0.87 


J98 


S03a 


0.012 


10.21 ± 0.37 


-11.91 ± 0.57 




7.74 ± 0.20 


P 


S08#l 














ESO 473-0024 


7.59 ± 0.06 


D 


S08 


c r\ A . 1 1 r» 

5.94 ± 1.19 


J98 


S03a 


0.011 


13.99 ± 0.81 


-14.85 ± 0.92 




7.45 + 0.20 


P 


S08 
















7.45 ± 0.03 


D 


S03b#2 














AM UlUo-ioz 


T T/i _i_ n on 


o 
r 


S08 


6.13 + 1.23 


TOQ 


oOJa 


0.00/ 


1 A _i_ f\ ^< 
14. 3o ± O.OJ 


1 A <0 _i_ n HQ 

— 14.j2 ± U. /o 




7.58 ± 0.04 


M 


L03(a) 
















7.62 ± 0.08 


P 


L03(a) 














NGC 625 


8.37 ± 0.06 


D 


S08#2 


3.89 ± 0.22 


RGBT 


C03 


0.010 


9.01 ± 0.26 


-18.95 ± 0.28 




7.69 ± 0.20 


P 


S08#2 
















5.19 ± 0.02 


D 


S03b#5 














ESO 245-G005 


8.17 ± 0.21 


D 


S08(c) 


4.43 ± 0.45 


RGBT 


K03 


0.009 


10.99 ± 0.50 


-17.24 ± 0.54 




7.61 ± 0.10 


D 


S08(d) 
















7.96 ± 0.20 


P 


S08 
















7.65- 8.20 


M 


M96, H03 




















M81 group 












DDO 42 


7.61 ± 0.06 


D 


S08 


3.20 ±0.41 


RGBT 


K02 


0.021 


10.13 ± 1.00 


-17.40 ± 1.03 




7.73 + 0.20 


P 


S08 
















7.9 ±0.01 


D 


197 














DDO 53 


7.83 ± 0.20 


P 


S08 


3.57 ± 0.25 


RGBT 


K02 


0.021 


12.64 ± 1.30 


-15.14 ± 1.31 




8.00 ± 0.23 




HH99(b) 














UGC 4483 


7.61 + 0.08 


D 


S08 


3.40 ± 0.20 


RGBT 


IT02 


0.020 


13.68 ±0.56 


-13.99 ±0.57 




7.43 ± 0.20 


P 


S08 
















7.56 + 0.03 


D 


V06 















Table 8. Abundances and NIR luminosities for the fiducial sample of galaxies. Distances are expressed in Mpc, and columns 6 
and 7 indicate the method used to compute the distance, and its reference. Our abundances have been corrected by -0.1 dex. For 
comparison we add also I Zw 18 (H = 16.03 ± 0.15; Thuan fT983]) and SBS 0335-052 (H = 16.09 ± 0.03, 15'' aperture; Vanzi et al. 
[20001) 
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luminosity enable a test of this method: they are UGCA 442, 
NGC 625, and ESQ 245-G005, and their distances obtained with 
the J98 method are respectively 18% smaller, 10% larger, and 
15% smaller than the ones quoted in Table|7l For the M81 group, 
all four galaxies have distances obtained with the RGBT tech- 
nique, and the DDO 42 distance has been measured by both 
Karachentsev et al. (l20Q2l K02) and Thuan & Izotov (l2Q05l 
TI05). Although the TI05 distance (3.42 ± 0.15 Mpc) has a 
smaller formal error, we adopt the K02 distance, as, doing so, 
we have consistent distances for three galaxies in our sample. 

All galaxies having direct abundance determinations are 
plotted in the upper panel of Fig. [S] A correction of -0.1 dex 
has been applied to our data, in order to be consistent with 
Skillman et al. (I2003bl) . There appears an obvious trend of in- 
creasing metallicity with H luminosity, except for the two galax- 
ies at 12 -h log(0/H) ^ 7.5 and Mr ~ -17 mag. These are 
ESQ 245-G005 and DDO 42: the l atter ha s an independent oxy- 
gen abundance from Izotov et al. (11997b , which is higher than 
our value and consistent with that of ESO 347-G017 at com- 
parable luminosity. Since the work of Roy et al. (119961) , it is 
known that DDO 42 displays a range of abundances consistent 
with a shallow gradient, so our low value probably reflects the 
external location of our target H ii region. In the following we 
will then adopt the Izotov et al. value which corresponds to a 
central region, and we note that a key point in the definition of 
an L-Z relation is that, for each galaxy, regions close to their 
centers should be selected. The case of ESO 245-G005 is simi- 
lar to that of DDO 42: in this galaxy, we measured both an ex- 
ternal and a central region, finding an abundance diff'erence of 
^0.5 dex. Our direct abundances (the first ones for this galaxy) 
confirm the metallicity gradient found by Miller (119961 M96) 
with an indirect method: our range is 12 -h log(0/H) =7.61 to 
8.17 while M96 finds a range 7.65 to 8.20. A similar result has 
been obtained by Hidalgo-Gamez et al. (I2003t . Again, we adopt 
the central abundance as the reference for this galaxy. 

From the sample plotted in the upper panel of Fig. O we 
extracted the galaxies listed in Table [H which are the fiducial 
objects defining our L-Z relation. Whenever possible we select 
abundances from our work, except for two galaxies: UGCA 442 
was not observed spectroscopically, and our observed DDO 42 
region was excluded for the reason discussed above. Our fiducial 
L-Z relation is plotted in the lower panel of Fig. [6l The obvious 
abundance trend with luminosity is confirmed by a correlation 
coeflftcient r = -0.93: a weighted least-squares fit yields the fol- 
lowing relation: 

12 -h log (0/H) = -0.20 ± 0.03 Mh + 4.54 ± 0.49 

with a dispersion of 0.1 1 dex around the fit. 

Table [8]includes also data for the blue-compact dwarf (BCD) 
galaxies I Zw 18 and SBS 0335-052. These are possibly the most 
metal-poor galaxies in the local Universe, so it is worthwhile to 
see how they compare with our sample of dwarf irregular galax- 
ies. This is done in Fig. [51 where the two objects are plotted 
with open symbols. They fall outside the trend defined by dirr 
galaxies, in the sense that their oxygen abundances are too low 
for their luminosities. There might be several explanations for 
this fact: the first possibility is that the gas enrichment of BCD 
galaxies is slower than that of dIrr galaxies, either because they 
produce less metals per stellar generation, or because they lose 
metals into the IGM more easily. So, as the luminosity of a BCD 
increases, its oxygen abundance does not increase as much as 
that of a dIrr of comparable mass. A second possibility is that 
we are not measuring the central abundances of the two BCDs. 
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Figure?. Comparison of the luminosities and abundances of 
our fiducial galaxies (filled diamonds) with the predictions of 
a set of closed-box models of increasing total masses (dotted 
line curves). For each model the evolution has been truncated 
when the remaining gas mass drops below 3x10^ At©, and the 
sequence of the end-points is outlined by the solid curve. The 
open symbols are the same as in Fig.[6j 



We have shown above that gradients as large as 0.5 dex are de- 
tected in dwarf galaxies, so the low abundance of the H ii regions 
in I Zw 18 and SBS 0335-052 might be explained by their non- 
central location. Indeed the H i maps presented in van Zee et al. 
(1998) and Hirashita et al. (2002 ), and compared to the optical 
and Hof maps, show that the star forming regions are slightly 
off'-center. However a counter-argument for I Zw 1 8 is that the 
galaxy is too small to have appreciable metallicity gradients. 
Finally in the case of SBS 0335-052, one might think that its 
distance is overestimated: if the distance were ~ 4 times smaller, 
then its luminosity would agree with the faintest objects in our 
sample. However at the relatively large recession velocity of the 
galaxy (~ 4000 km sec"^) peculiar motions cannot induce such a 
factor of four uncertainty in the distance (e.g., Tonrv l2000l) . The 
conclusion is then that these two BCD galaxies must have had 
diff'erent evolutionary paths than dIrr galaxies. 

This remains true even when the two galaxies are compared 
to the general population of BCDs. Shi et al. (2005 ) have shown 
that the L-Z relation of BCD galaxies is similar to that of dIrr 
galaxies, so at the luminosity of I Zw 18 and SBS 0335-052 a 
typical BCD has higher oxygen abundances than these two ob- 
jects. On the other hand the Shi et al. L-Z relation shows a very 
large scatter: this is not quantified by the authors, but in their 
Fig. 3 one can see object-to-object abundance diff'erences of up 
to 1 dex. Shi et al. propose that the scatter might be due to dif- 
ferences in SFH or the evolutionary status of the current star- 
burst, or to diff'erent starburst-driven outflows or gas infall rates. 
The abundances of I Zw 18 and SBS 0335-052 are then partic- 
ularly low, but still compatible with those of BCDs in general, 
and could be due to one of these processes. 

5. Discussion 

5.1. Origin of the lumlnoslty-metaHlclty relation 

As recalled in the introduction, in the case of dE/dSph galax- 
ies the L-Z relation is interpreted in the scenario of mass-loss 
through galactic winds. The gas escapes more easily from low 
mass galaxies, so their chemical evolution is truncated before 
that of more massive galaxies. Some previous studies have pro- 



posed that dSphs and dirrs share the same L-Z relation, thus 
extending the mass -loss s cenario to irregulars (e.g. Skillman, 
Kennicutt, & Hodge [T989l) . However, recently it has been shown 
that the abundances of dIrrs are lower than those of dSphs 
of comparable luminosity, whether we consider their gas con- 
tent (e.g. Mateo 1998) or their RGB stars (e.g. Grebel 2004). 
Therefore, if we want to retain the mass-loss scenario for both 
classes of dwarfs, then we must assume not only that this pro- 
cess leaves behind an L-Z relation after the evolution has stopped 
(dSph), but also that it induces a (possibly different) relation very 
early in the history of galaxies (dirr), and is able to maintain it 
throughout their lives. 

More input on the issue of the chemical evolution of dIrrs has 
been provided recently by Skillman et al. (I2003bt and Pilyugin 
et al. (12004b . Both studies conclude that the chemical evolution 
of dIrrs can be approximated by a closed-box, provided that 
an effective oxygen yield ^ 1/3 of the theoretical value (e.g., 
Maeder 119921) is adopted. A low yield can be the signature of 
gas exchange with the environment, and this is indeed observed 
in a few cases, as discussed in the introduction. If this conclusion 
is valid, then it is hard to imagine how mass-loss can produce 
the L-Z relation, because the galaxy mass is not a parameter of 
the closed-box model. Conversely, one would expect that the ef- 
fective yield does depend on the mass of the galaxy, since the 
efficiency of mass-loss via galactic winds should depend on the 
depth of the potential well. Note also that Lee et al. (2006 ) have 
found a few dwarf galaxies with large yields, a result difficult to 
reconcile with the galactic wind hypothesis. In conclusion, there 
does not appear to be any obvious interpretation of the observa- 
tional facts and the situation remains confusing. 

Yet, pushing further the hypothesis of a low and univer- 
sal yield, we can gain some insight by looking at the dotted 
curves in Fig. (T] which are in fact closed-box models com- 
puted with constant yield P(o/n) = 1.6 x 10"^, and total masses 
mtot = ^gas +^stars Varying between 4x10'^ At© and 1 .4x 10^ At©. 
We also assumed M/Lh = 1.2 (typical of a globular cluster-like 
population, e.g. Bruzual & Chariot 120031) . and Mh,q = 3.1 mag 
(Astrophysical Quantities); the only parameter varying along a 
track is the gas mass fraction ji = mgas/mtot- If one assumes that 
all galaxies were born at the same time, then this plot is telling us 
that more massive galaxies evolve faster along their tracks, i.e. 
dji/dt is larger for more massive galaxies. If the galaxies were 
formed at different times, then larger galaxies had more time to 
build up their metals. Skillman et al. (12003 ab computed SFRs for 
the Sculptor group galaxies: if their data are plotted in the form 
djd/dt = SFR/(mgas + mstars) vs. (mgas + mstars), then a general 
trend of d/i/dt increasing with total mass can be seen. This con- 
firms that larger galaxies are more effective in converting their 
gas into stars. 

The conclusion is that, while a galaxy is still evolving (such 
as a dIrr) the L-Z relation could be a combined effect of mass- 
loss and of more efficient gas processing. A similar conclusion 
was reached by Pilyugin & Ferrini (12000b : according to them, 
the L-Z relation is a combined effect of smaller gas-loss and 
higher astration level as the mass increases. Another prediction 
of this scenario is that the slope of the L-Z relation should in- 
crease with time, and the effect would be even larger if the gas 
is removed at earlier stages in smaller galaxies. The reason why 
astration becomes more efficient as the galaxy mass increases, 
remains to be understood. On the other hand, we should recall 
that one of the most accepted facts ab out sta r formation is that 
it depends on the gas density (Schmidt l 1959b , a conclusion also 
suggested by recent simulations (e.g., Chiosi & Carraro [2002b . 



So it might be that as its mass increases, a galaxy becomes more 
effective at compressing its gas. 

The open symbols in Fig. [7] represent again the two BCD 
galaxies I Zw 18 and SBS 0335-052. In Sect. |4] we concluded 
that they do not follow the trend defined by dIrr galaxies because 
their chemical evolution must be different. Our simplified ap- 
proach would tell us that I Zw 1 8 is on the track of a ~ 3 x 10^ At© 
galaxy, while SBS 0335-052 is on that of a ~ 2 x 10^ At© object, 
and that - compared to a dIrr of similar total mass - they did 
not convert much of their gas into stars, thus having very low 
metallicities. However it is also possible that their evolutionary 
paths are different: for example their compactness might imply 
higher densities, hence SFRs, which in turn would trigger more 
powerful stellar winds and remove metals more easily than in 
dIrr galaxies. Or they might be unusually young. 

The closed-box scenario itself does not make any prediction 
about the final abundance a galaxy can reach: actually it predicts 
arbitrarily large values of Z as the gas mass fraction tends to 
zero, since Z = p ln(l/yu) (Searle & Sargent 1972). However, if 
we assume that the mass of each new generation of stars, msF, is 
roughly constant, then the maximum metallicity we can measure 
will be Zmax = P ln(mtot/^SF) ^ P In mtot, i-e. it will be propor- 
tional to the mass of the galaxy. This is shown in Fig. [71 where 
we assumed msF = 3 x 10^ At©, a mass of stars that can be created 
in ~ 10^ yr for a typical SFR -0.1 At©yr"^ This formulation 
seems to provide a natural upper limit to the measurable metal- 
licities, and it works even for large galaxy masses. Indeed, the 
dotted line in the top panel of Fig. [5] represents the same rela- 
tion extended to mtot = 2 x 10^^ At©, and with the exception of 
a couple of objects, the metallicities of star-forming galaxies are 
always lower than this limit. It is quite surprising that an upper 
limit to the metallicity is provided by such a simple expression, 
as it is an oversimplification of galactic chemical evolution. Still, 
it might suggest that a characteristic mass is involved in the star 
formation process. 

5.2. Comparison with other galaxies 

The literature offers only a few 0/H measurements for galaxies 
as faint as those in our sample, so a direct consistency check of 
our results is not possible. However, one can find a few dIrrs 
with oxygen abundances and total infrared magnitudes allowing 
us to verify our relation statistically. 

In order to minimize any systematic errors, we limited our- 
selves to using the infrared magnitudes from three sources: the 
2M ASS E xtended Source Catalog (Ja rrett e t al. 12000b . Cairos et 
al. (12003b . and Hunter & Elmegreen (12006b . No corrections for 
the different filter systems were applied. Furthermore, we used 
only oxygen abundances derived with the temperature sensitive 
method. The only exception is DDO 43, where T^^ 10000 K was 
assumed instead of being measured. Finally, the distances for 
galaxies closer than ~10 Mpc were based on direct measure- 
ments rather than on the Hubble flow. The data used for the 
comparison are listed in Table [51 More than one 0/H measure- 
ment per galaxy is listed, whenever available. Radial velocities 
are taken from the NASA extragalactic database. If the literature 
sources do not list the 0/H uncertainty, 0.1 dex was adopted. 
The distances flagged with "h" are based on the Hubble flow 
(Ho =72 km s"^ Mpc"^). The //-band extinction Ah has been 
taken from Schlegel et al. (1998). The table data are presented 
in graphical form in Fig. [8] and compared to our L-Z relation. It 
appears that the properties of our targets, dIrrs in the Sculptor 
and M81 groups, are generally similar to those of nearby dIrr 
galaxies, but the L-Z relations yielded by the two sets of data are 
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Figure 8. Abundance vs. absolute H luminosity for galaxies with 
data taken from the literature (filled triangles) and for our tar- 
gets (filled circles). In the upper panel, the dashed line is an un- 
weighted fit to the literature data, while our fiducial relation is 
shown in the lower panel. In both panels dotted lines outline the 
r.m.s. dispersion around the fit. 



radically diff'erent. The slope of an L-Z relation based on litera- 
ture data is 2.5 times smaller than the one we obtain (-0.08+0.02 
vs. -0.20 ± 0.03), and its dispersion is 0.18 dex, i.e. 58% larger 
than the scatter around our relation (0.1 1 dex). 

The large scatter in the literature data might be ascribed to 
the uncertain distances, and it is also likely related to the ran- 
dom location of the H ii regions considered inside each galaxy. 
Moreover, this large scatter seems to aff'ect faint galaxies more 
than luminous ones, since all objects brighter than Mh - -18 
exhibit a smaller scatter. Note that the literature compilation in- 
cludes galaxies in diverse environments, another factor which 
might account for the larger spread in the literature L-Z relation 
in comparison with our L-Z relation for the group galaxies. 

Finally, although we tried to use a sample that is as homoge- 
neous as possible, and the method for the abundance calculation 
is formally the same in all these works, yet systematic diff'er- 
ences are still possible (as demonstrated in Sect. 13.11) . and they 
can explain part of the scatter. A detailed discussion of all the 
single galaxies contained in Table [9] is beyond the scope of this 
paper, but the conclusion is clear: any quantitative interpretation 
of the L-Z relation based on literature data should be performed 
with great caution. However, this conclusion might be revised as 
more and more large and homogeneous data samples appear in 
the literature. 

5.3. A dwarf vs. giant galaxy dichotomy? 

The unique study of the L-Z relation for giant galaxies that in- 
cludes NIR, //-band magnitudes is that by Salzer et al. (120051 
hereafter S05). They obtained spectra of emission-line galax- 
ies in fixed apertures of V!5 or T\ and computed oxygen abun- 
dances with a reduced number of emission lines, because of the 



limited spectral coverage of their data: using additional spectra, 
they first derived metallicities for a subsample of galaxies using 
both the direct and the empirical methods, and then calibrated 
them against [Nii]/16583/Ha and [O iii]^5007/HyS. The empir- 
ical method they used is the Pilyugin (2000; POO) calibration 
for the lower branch of the 12 + log(0/H) vs. 7?23 relation, and 
three calibrations for the upper branc h: Edm unds & Paget (119841 
EP), Kennicutt, Bresolin, & Garnett (l2003l KBG), and Tremonti 
et al. (12004b . In Fig. [9] we have plotted our data together with 
Salzer's et al. data. The figure reveals that the scatter in the L-Z 
relation for giant galaxies is very large, compared to that of our 
dirrs, perhaps due to fixed-aperture eff'ects and uncertainties spe- 
cific to the empirical methods (see also the discussion in S05). 
The straight lines in the upper panel are the fits obtained by S05 
and by us: taken at face value, this plot would suggest a well- 
defined off'set between the L-Z relation of giants and dwarfs. The 
slope of the S05 L-Z relation is in fact close to what we find for 
dwarfs alone, namely -0.215 ± 0.003 and -0.201 ± 0.004 for 
abundances obtained with the EP and KBG calibrations, respec- 
tively. Finally, to be consistent with S05, in the central panel we 
have plotted our abundances computed with the same POO cal- 
ibration of the empirical method, allowing two galaxies to be 
added (see abundances identified by "P" in Table [7]). The off'set 
between galaxies in the two mass ranges seems to disappear: it 
is replaced by what looks like diff'erent L-Z relations for dwarf 
and giant galaxies, with a separation at Mh ~ -20. A fit to S05 
galaxies fainter than this limit yields the dotted line shown in the 
panel, while the solid line is a fit to our dwarfs. The two relations 
have a similar slope, and a ~ 0.3 dex off'set, with the dwarfs of 
our sample being on average more metal poor than those of S05. 
This might result from the two-step calibration of S05, which 
introduces additional uncertainties in the empirical method - a 
method which, in any case, does not guarantee abundance preci- 
sions better than 0.2 dex -. The most important thing to notice, 
however, is that extremely diff'erent L-Z relations are obtained 
for the dwarfs of our sample, if one uses abundances computed 
by the two methods, direct or empirical. The solid line in the up- 
per panel of Fig. [9] (corresponding to direct abundances in the 
region of dwarf galaxies) has a slope of -0.05 dex mag"^ and an 
r.m.s. dispersion around this fit of 0.14 dex. The scatter is larger, 
and the slope is four times smaller when empirical abundances 
are used for the same sample of dwarfs (central panel). Striving 
for abundances obtained with the direct method requires a much 
larger observational eff'ort, but the eff'ort is wholly justified by 
the better constraints it provides on galaxy evolution. 

Interestingly, using direct abundances from the literature 
yields an L-Z relation similar to that obtained by fitting the 
S05 dwarfs. This can be seen by comparing the solid line 
in the bottom panel of Fig. [51 to the dotted line in the cen- 
tral panel. The zero points are similar, and the two slopes are 
-0.08 ± 0.02 dex mag"^ and -0.04 ± 0.02 dex mag"\ respec- 
tively. The conclusion is that the use of direct abundances is not 
a sufficient condition: distances must be know with great accu- 
racy as well. In that respect, let us point out that selecting a sam- 
ple of dirr galaxies in groups has been a winning strategy. Even 
assuming just an average distance for all galaxies in a group has 
allowed us to obtain essentially the same L-Z relation (Saviane 
etal. [20041 [20051) . 

The question of a dwarf-giant dichotomy is then open: mov- 
ing from the empirical to the direct abundances, the L-Z rela- 
tion for dwarfs becomes much better defined, parallel to the one 
obtained for more massive galaxies, and with a substantial off'- 
set. It remains to be seen what would happen to the L-Z rela- 
tion of giant galaxies if one uses direct abundances: an impor- 
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Figure 9. Oxygen abundances vs. H luminosity for the Salzer et al. ('2005') emission-line galaxies (KBG calibration), our dirrs, 
and data compiled from the literature. In all panels the dashed line is the fit obtained by Salzer et al. to their data points (small 
open circles). In the top panel, the thin solid line is the same curve as that shown in Fig. [71 while the filled triangles and solid line 
correspond to our fiducial L-Z relation. In the central panel, the dotted line is a weighted fit to Salzer et al. galaxies fainter than 
Mh = -20, while the filled triangles and solid line relate to our L-Z relation obtained using abundances re-computed by us with the 
empirical, indirect method. In the bottom panel, we show the L-Z relation from literature data (filled triangles and solid line). 



tant extension of this project would be to measure direct abun- 
dances for at least some of the S05 galaxies, and check whether 
the observed off'set is confirmed. Since the SOS sample includes 
star-forming galaxies, one might suspect that some of the giant 
galaxies are shifted in lumino sity b ecause of the presence of a 
star-burst: however, Lee et al. (120041) have found that this shift is 
a few tenths of magnitude in the 5-band, compared to quiescent 
galaxies (however see also Bica et al. 1990). We expect that the 
eff'ect in the IR must be even lower, and certainly not comparable 
to the ~ 2 magnitudes we observe. 

A dwarf-giant galaxy dichotomy has been discussed, e.g., 
in Garnett (120021 hereafter G02), where the trend of eff'ective 
yields = Z/lnQi'^) vs. galaxy luminosity or rotational ve- 
locity is studied. It is shown that dwarf galaxies (those with 
Vrot ^100 km s"^) have lower eff'ective yields than massive spiral 
galaxies, which is interpreted as due to gas outflows in dwarfs. A 
similar conclusion has been reached by Tremonti et al. (120041) . 



In G02 the stellar mass that enters the gas mass fraction ji was 
computed from blue luminosities, assuming a color-M/L rela- 
tion. Since more accurate stellar masses can be computed with 
NIR luminosities, the availability for massive spiral galaxies of 
abundances obtained with temperature- sensitive methods would 
allow more accurate efl'ective yields to be computed as well, and 
therefore to update the G02 study. In particular it would be in- 
teresting to check whether indeed dwarf galaxies can lose up to 
90% of their metals and be the main source of enrichment of the 
inter-galactic medium. 



6. Summary and conclusions 

We have measured oxygen abundances in H ii regions located in- 
side a number of dirr galaxies, belonging to the nearby Sculptor 
and M81 groups. The abundances were measured with the 
temperature- sensitive (direct) method, based on the ratio of the 



auroral line [O iii].14363 to the nebular lines [O iii].1.14959, 5007, 
and were complemented by direct abundances from the literature 
for two additional galaxies. The weak [Oiii]/14363 line could 
be measured in Hii regions belonging to five galaxies of the 
Sculptor group and two galaxies of the M81 group. Metallicity 
gradients were detected in ESO 245-G005 and DDO 42, con- 
firming earlier findings, so only the central highest-metallicity 
regions were considered for the derivation of the luminosity- 
metallicity relation. This forced us to discard our DDO 42 mea- 
surement, and adopt a literature value. Our fiducial sample then 
includes six dirr galaxies of the Sculptor group and two galaxies 
of the M81 group. 

For these eight galaxies we have obtained deep NIR, if -band, 
imaging which allowed us to perform surface photometry and 
compute their total luminosities. Thanks to the availability of 
distances with < 15% errors, the galaxies could be placed in 
the 12 -h log (0/H) vs. Mh diagram, revealing a clear L-Z rela- 
tion with a small 0.1 1 dex scatter around the average trend. The 
scatter is smaller than that of relations obtained at o ptical w ave- 
lengths (e.g. 0.161 dex using 5-band data, Lee et al. 12006b . and 
is comparable to the one obtained by Lee et al. (12006b with their 
MIR [4.5yum]-band data (0.12 dex). Assuming the existence of 
a fundamental mass-metallicity relation, the improved definition 
of the L-Z relation at NIR luminosities must be due to the fact 
that it better traces the underlying relation with mass, since the 
NIR mass-to-light ratio is more sensitive to the integrated star 
formation history of a galaxy. On the contrary, the blue mass-to- 
light ratio is more sensitive to the instantaneous SFR, which has 
a large galaxy-to-galaxy scatter (Saviane et al. 120041 Salzer et al. 
[20051) . 

Our work and that of Lee et al. (12006b are the ones that man- 
aged to obtain the best defined L-Z relation, and as an additional 
advantage, our standard NIR band allows an easy comparison to 
other relations obtained in independent studies. Indeed we com- 
pared our relation to that of Salzer et al. (12005b . who are the only 
authors providing an //-band relation, at the same time extend- 
ing it to giant galaxies. Unfortunately a direct comparison with 
our accurate abundances is not possible, since the abundances 
of S05 are obtained with the so-called strong-line method (indi- 
rect or empirical method). To be consistent with S05 the abun- 
dances of our dwarf sample were also computed with the empir- 
ical method, and although the new L-Z relation has more scat- 
ter than the one using direct abundances, in this way our dwarf 
galaxies can be placed in a same graph together with the galaxies 
analyzed by S05. The slope and zero-point of the relation for our 
dwarfs, using indirect abundances, are diff'erent than those of the 
relation obtained with direct abundances, and, when compared 
to S05 giants, a possible break-point appears at Mh ~ -20. 
While this could suggest a diff'erent gas-consumption mecha- 
nism for dwarf and giant galaxies, a solid conclusion cannot 
be established at this stage. Indeed, the version of the empiri- 
cal method employed by S05 is diff'erent from that classically 
used, and so a slight inconsistency between our abundances and 
those of S05 may still exist. The best way for giving the final 
word on this possible giant-dwarf dichotomy would be to obtain 
direct abundances for a good number of the KISS galaxies. And 
precisely, such a project has been started recently by some of us 
(see Saviane et al. 12007b . 

One of the motivations for the current study was to rem- 
edy the lack of homogeneous abundances and luminosities in 
the literature, but still it was interesting to make the experiment 
of assembling an L-Z relation based on literature data, and see 
whether our approach really turned out to be superior. We con- 
cluded that indeed mixing data from a variety of sources builds 



an L-Z relation aff'ected by large uncertainties, and very diff'erent 
from the one obtained with a controlled sample and controlled 
methods. Inferences based on literature data are thus to be taken 
as qualitative at best, and should be discarded, if possible. 

Limiting ourselves to our sample of dirrs, we attempted to 
explain the L-Z relation assuming that the chemical evolution 
of these objects is similar to that of a closed-box model with 
an eff'ective yield which is 1/3 its value in the solar vicinity 
(Skillman et al. I2003bl Pilyugin et al. l2004l) . If this approxi- 
mation holds true, and if we assume that the chemical evolu- 
tion has started at the same time irrespective of mass, then the 
conclusion is that more massive galaxies have faster chemical 
evolutions (gas-consumption rates). Alternatively, more massive 
galaxies had more time to build up their metal content. 

Due to a number of reasons (mainly adverse weather condi- 
tions) we could assemble a large database only for the Sculptor 
group of galaxies, while only two galaxies of the M81 group en- 
ter our L-Z relation. Although these two dwarfs seem to follow 
the same relation as that defined by the Sculptor group galaxies, 
more M8 1 dIrrs need to be measured, before one can understand 
whether the higher density in that group is influencing the chem- 
ical evolution of its members. 

Indeed in general terms our L-Z relation is based on a small 
number of galaxies, therefore more objects need to be added to 
Fig.[6]in order to confirm our L-Z relation with a larger sample. 
For example looking at the number of dIrr presented in Lee et al. 
(12006b . we might be able to add ~ 20 galaxies to our database. 
And lOm-class facilities will have to be used to reach this goal 
in a reasonable time. In the short term we plan to keep collect- 
ing data for the M81 group, and start a similar study for the 
Centaurus A group. A natural extension of the project would be 
to investigate clusters of galaxies, but this is a difficult task from 
the observational point of view, since the two nearest galaxy 
clusters are Virgo and Fornax, at 22 and 24 Mpc respectively 
(Ferguson & Sandage [19901 hereafter FS90). A more viable al- 
ternative would be to observe the Leo group of galaxies. At a 
distance of 18.2 Mpc (FS90), this group is the fourth nearest one, 
and has a relatively high velocity dispersion of 250 km sec~^ . Its 
properties (density, number of galaxies, dwarf-to-giant galaxy 
ratio, etc.) are intermediate between those of nearby groups and 
those of clusters of galaxies (see e.g. Figure 1 and following in 
Ferguson & Sandage 1991). It would then allow us to begin the 
exploration of a cluster-like environment, with NIR cameras at 
lOm-class telescopes. 
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Appendix A: Reduction of EF0SC2 long-slit 
spectroscopy 

AA. Extraction of spectra 

All nights were treated separately, and reductions were carried 
out using a customized version of the EF0SC2 quick-look tool 
(http : //www . Is . eso . org/), which is based on the eso-midas 
data reduction system. The main reduction steps are summarized 
as follows. 

All bias frames were inspected to ensure that their level was 
the same, and an average bias was then computed and subtracted 
from the rest of the frames. The master flatfield was computed 
by taking the median of five normalized frames and then sub- 
tracting the spectral signature of the quartz lamp. The latter was 
computed by first averaging the flux along the spatial direction, 
then fitting a polynomial of degree 6 along the dispersion, and 
finally expanding it back into a bidimensional frame. The wave- 
length calibration was computed based on the spectrum of the 
internal Helium- Argon lamp, taken with the wide slit. Since 
the projection of the 5'' wide slit on the CCD is in the same po- 
sition, the calibration was used both for the spectra of the H ii 
regions and the spectrophotometric standards. A transformation 
of the form (x, y) (A, s) was found, thus correcting also for the 
distortion of the spectra introduced by the optical system. The 
transformation was modeled with polynomials of degree 3 and 2 
along the dispersion and the spatial direction, respectively. After 
applying the calibration, the spectra were linearly rebinned with 
a constant step of AA = 4. 23 A. The bidimensional sky frame 
was created by first sampling the sky spectrum in two windows 
flanking the target spectrum, and then fitting the spatial gradient 
with a polynomial of degree 1 . After subtracting the sky frame, 
the extraction window of the emission-line spectra was decided 
by looking at the spatial profile of the Ha line and making some 
experiments until the best S/N was found. On the other hand, 
the flux of the spectrophotometric standard stars was summed 
over almost the entire spatial profile, leaving out just the two sky 
windows. As it was recalled above, the response function for 
the flux calibration was computed using the star LTT 3218. To 



do this, the instrumental spectrum was corrected for atmospheric 
extinction, divided by the exposure time, and then divided by the 
pubHshed spectrum, and the ratio was fitted with a polynomial 
of degree 12. The spectrophotometric standard was observed at 
several airmasses, and we checked that the extinction-corrected 
and normalized spectra are almost identical, thus confirming the 
extinction curve for La Silla provided by midas. The response 
function was then computed as the median of the five func- 
tions found for the two nights. The pre-processed, wavelength- 
calibrated, sky- subtracted and extracted spectra were thus finally 
flux-calibrated using this response function. 

A.2. Flux measurements 

The fluxes were measured within the alice contex in midas. 
A fourth-degree continuum was fitted to the line-free spec- 
tra, and then for each line the flux was summed within two 
wavelengths bracketing the line. The blends H7-h[Oiii]/14363, 
HQf-F[Nii]^6583 and the [Sii] doublet ^^6716,6731 were de- 
convolved using Gaussian fits, and we ensured that the sum of 
the single line fluxes equaled that of their blend. The measured 
line ratios were corrected for the efl'ects of reddening, measured 
by comparing the observed and the expected hydrogen l ine ra- 
tios. The theoretical H line ratios of Hummer & Storey (119871) 
were used, assuming an electron temperature Tq ^ 10000 K. 
We set R = hUnp, Ro = /i,o/%,o and if(A) = [f(A) - f(Hp)] 
so that C = [logRo/R]/(f(A). The reddening constant was then 
computed using the Cardelli et al. (119891) extinction law, where 
f(A) =< A(A)/A(V) >, and we assumed Ry = 3.1. We obtained 
the value of C both from the ratio Ha/HyS and from HjS/Hy. These 
line ratios are afl'ected by the unknown stellar absorption under 
the H emission lines: we then corrected the fluxes assuming a 
constant EW of the absorption lines, and adopted the EW that 
yielded the closest values for the two C. In some cases we were 
not able to resolve Ha from the [Nii] lines; however, the flux 
of the latter lines is at most a few percent of that of Ha, so we 
assumed that this correction is negligible. Since the Ha/UjS ra- 
tio is the more accurate, the value from this ratio was adopted 
to correct the observed fluxes, which were finally normalized to 
UJ3 whose flux was set to 1 . 



Appendix B: Reduction of KAST long-slit 
spectroscopy 

B.1. Extraction of spectra 

In the following we will call blue, red, full, and high-resolution, 
the spectra taken with the blue arm and grism 830/3460, with 
the red arm and grating 300/4230 (plus dichroic D46), with 
the red arm and no dichroic, and with the red arm and grating 
1200/5000, respectively. 

The frame pre-processing, wavelength calibration, and flux 
calibration was done within eso-midas (and its "context" long), 
with a modified version of the quick-look tool described above. 
The procedure was similar to the EF0SC2 data reduction, ex- 
cept for the bias subtraction, which is done directly by the Kast 
acquisition system. All the frames were trimmed to remove the 
overscan regions, and then master flat-field (FF) and arc images 
were created. Three to five dome flat-field frames per spectro- 
graph configuration were available, so a median frame was first 
created, and then the spectral signature of the lamp was removed 
by fitting a polynomial surface to the median FF, assuming that 
the polynomial does not vary along the spatial direction. 
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Figure B.l. Response functions for each date of the Lick run: 
the solid curves are the functions for the blue and red spectra, 
the dashed curve is the full spectrum function (no dichroic), and 
the dotted lines are the functions for the high resolution spectra. 
The functions have been normalized to the peak response. 

Table B.l. Dispersions and wavelength range for each date and 
Kast configuration 



Kast setup disp. wavelength range 

Arm/dichroic/gratinglgrism [A/px] 





Dec 17, 2001 






blue/D46/830 3460 


1.1 


3198.0 


4532.2 


red/D46/1200_5000 


1.2 


5856.5 


7240.8 


red/D46/300_4230 


4.6 


4216.2 


9099.2 


red/none/300_4230 


4.6 


3688.1 


9097.7 




Jan 14, 2002 






blue/D46/830 3460 


1.1 


3213.1 


4552.3 


red/D46/1200_5000 


1.1 


5874.4 


7194.9 


red/D46/300_4230 


4.6 


4170.7 


9052.5 




Jan 15 2002 






blue/D46/830_3460 


1.1 


3213.3 


4552.7 


red/D46/300_4230 


4.6 


4149.0 


9032.9 



Our main concern in the data reduction was to ensure that, 
in the overlap region, the fluxes of blue and red spectra are cal- 
ibrated with the smallest systematic uncertainties. This was par- 
ticularly difficult for the shortest wavelengths of the red spec- 
tra (blueward of ~ 5000 A), where the response function drops 
steeply from 50% to zero in just ~ 300 A (see below and 
Fig. IB. II) . This gradient is extremely sensitive to the FF correc- 
tion in this spectral region, and it was found that a polynomial of 
degree 7 gave the best response function for the red spectra. For 
all other spectra a degree 6 polynomial gave the best results. 

The arcs were used to remap the bidimensional spectra from 
the pixel-pixel space to the wavelength-pixel space, thereby also 
correcting the off'-axis distortions of the optics, such that the 
spectrum could be extracted by just defining a window in the 
spatial direction. The two Kast lamps are filled with NeAr and 
HeHgCd, and since no cadmium or mercury line lists are avail- 



able in the midas system, a suitable line list was built from var- 
ious sources. The computed dispersions and wave leng th cover- 
ages for each Kast configuration are given in Table IB. ll The dis- 
persion is ~ lA px"^ for the blue and high resolution spectra, 
and 4.6A px"^ for the red and full range spectra. The resolution, 
as measured on the arc lines, is ~ 3 A FWHM and ~ lOA FWHM 
for the high and low resolution spectra. 

In the next step, all the spectra were divided by the master 
FF and calibrated in wavelength, and then a preliminary extrac- 
tion with a large window was done, in order to create sky images 
for each spectrum, that were used for the correction of the spec- 
trograph flexures (see below). The sky frames were constructed 
by first creating the sky spectrum on the two sides of the cen- 
tral object, as an average within two windows. The flanking sky 
spectra were used to fit the sky gradient in the spatial direction, 
assumed to be linear and independent of wavelength, and the 
bidimensional artificial sky frames were made by simulating this 
gradient. 

By comparing the position of sky lines in diff'erent spectra, it 
was soon realized that the spectrograph is subject to large flex- 
ures, with spectrum-to- spectrum difl'erences of up to 10 px, i.e. 
more than 40 A for the red spectra. This fact has several impli- 
cations, and at this point of the reduction it afl'ects the defini- 
tion of the response function, in particular at the shortest wave- 
lengths of the red spectra. For this reason, the standard star spec- 
tra were registered to rest- frame wavelengths before creating the 
response function. Due to the short exposure times, the sky lines 
are almost absent from the standard star spectra, so the shift in 
wavelength was computed using the central wavelength of Ha 
for the red spectra, and Hy for the blue spectra. This approach 
works if the eff'ect of radial velocities on the position of the stel- 
lar lines is negligible. This assumption was checked by compar- 
ing the relative position of Ha and the telluric triplet at ~ 6870A. 
The average diff'erence over all standard stars is A/1 = 309.1+0.7 
A, and since the scatter is much less than the spectral dispersion, 
we conclude that indeed the eff'ect of radial velocities (both in- 
trinsic and due to the solar motion) is negligible, and that the 
barycenter of the telluric triplet is at 6871.9 A. 

The shift in wavelength was then subtracted from the start- 
ing wavelength of the bidimensional spectra, the standard star 
spectra were extracted again, and the response function was de- 
fined by comparison to the published spectrum. Again, the most 
difficult task is following the response function for the red spec- 
tra blueward of ~ 5000A. It was found that a polynomial fit is 
not adequate even using the highest degree implemented in mi- 
das (30). We then used a spline of degree 1, which gave good 
results since the sample rate of the published fluxes is a few A. 
The master response function was decided by flux-calibrating 
all the standard stars with a given response function, and then 
by taking the one that gave the best average calibration for all 
standards and over a range of airmasses (see Fig. IB.lt . Due to 
the large wavelength coverage, the most critical case is that of 
the red spectra, where we expect a maximum difl'erence of ~ 5% 
from the red to the blue part of the spectrum. For these spec- 
tra, the response function is able to correct the instrumental flux 
down to 4300A. 

Finally the galaxy spectra were extracted and flux-calibrated 
with a procedure resembling that adopted for the standard stars, 
but in this case the flexure correction was done with the very 
strong sodium sky line at 5 892 A. In order to have a good sub- 
traction of the sky lines, the sky windows were chosen close to 
the object spectrum, and the best windows were decided on an 
object-by-object basis with several experiments. Also the win- 



dows for the object spectra needed a careful selection, in order 
to enhance the faint lines over the background, and to avoid as 
much as possible the places where the sky line subtraction was 
less than optimal. 

B.2. Creation of merged spectra 

Using the continuum in the overlap regions, the final spectra 
have been obtained by merging the blue, red (or full in the case 
of DDO 42), and high resolution spectra. This implies that a re- 
liable merging could be produced only for those spectra where 
a meaningful continuum could be defined (i.e. those with the 
higher S /N), namely for the regions in DDO 42, DDO 53, and 
UGC 4483. First, the continua have been equalized in the over- 
lap regions, and then the spectra have been rebinned to a com- 
mon lA step, which is smaller than any of the other disper- 
sions. The equalization of the continua has been done in the 
log(A) - log(flux) plane: a straight line has been fitted to the con- 
tinua, and then the diff'erence d of the two fits at the midpoint 
of the overlap region has been taken. The normalization factor 
is then 10*^. In the merged spectrum Ha belongs to the high- 
resolution spectrum, HjS belongs to the red (or full) spectrum, 
and Hy belongs to the blue spectrum. These hydrogen lines were 
then used to perform some health checks on the merged spectra. 
Both Ha and HJS are present in the red (or full) spectrum, so we 
first checked that the ratio of these two lines is the same in the 
merged spectrum as in the red (or full) spectrum. If that was not 
the case, corrections to the equalization factors were introduced. 
Hy is present both in the blue and in the red (or full) spectrum, 
however in that spectral region the response function of the red 
(or full) spectrum reaches almost zero, so its flux is rather un- 
certain. Comparing the flux of this line in the two spectra is then 
not really meaningful, so we followed another approach. For the 
Sculptor group dwarfs, Table|4]shows that the average flux of Hy 
is 0.46 X HJS, with a dispersion of 0.04 x H/S. Thus we made sure 
that no strong deviations from this value were obtained: if the 
Hy flux was more than 3cr discrepant from the quoted value, we 
adjusted the equalization factor to bring the flux into agreement 
with the average value while maintaining a good overlap of the 
continua. This meant reducing the deviations to less than 2cr. 

B.3. Flux measurements 

Line fluxes were measured in the way described in Sect. IA.21 
regarding EF0SC2 measurements. The difl'erence is that a dif- 
ferent set of lines had to be deblended. In particular, [O ii] 3727 
was deblended from two nearby unknown lines at 3722A and 
3736A, [O iii]4363 was deblended from Hy4340, and finally the 
Gaussian deblending was employed to separate [Oii]7320 from 
[Oii]7330. 

Appendix C: Reduction of near-IR imaging 

The first step of the data reduction was to remove the sky back- 
ground, and simultaneously, the dark current and bias. We sub- 
tracted from each object image the average of the preceding and 
the succeeding sky images. This strategy accounts for sky back- 
ground variations with characteristic time comparable or a few 
times longer than the time spent at each pointing. For the first 
and the last object, we could only subtract the nearest sky image. 
Next, we carried out the flat field correction. A common median 
zero point was imposed to all available sky images, a median 
was obtained to produce a flat field, which was normalized to 



unity. Then, all sky- subtracted object images were divided by 
the flat-field. The bad pixels were masked out and all object im- 
ages were aligned and combined to produce a final image for 
each target. In addition, for the SOFI data, we carried out an il- 
lumination correction to account for the diff'erence between the 
sky illumination and the illumination of the dome screen used 
for flat-fielding, as described on the instrument webpages. The 
entire data reduction was carried out with standard IRAF tasks. 

The photometric calibration was carried out by observations 
of standard stars from Hunt et al. (1998) for the Northern ob- 
jects, and from Persson et al. (11998) for the Southern ones. The 
INGRID data were photometrically calibrated with measure- 
ments of 129 standard stars in 13 fields from the night of Jan 
16, 2003, and the SOFI data were calibrated with 13 standard 
stars observed during all nights. Some galaxies were observed 
twice at the NTT to verify the self consistency of the photome- 
try and because the night of Oct. 15, 2002 was non-photometric, 
so all targets had to be re-observed during the next night with 
shorter integration times. Finally, the calibration of ES0473- 
024 is based on a single standard, observed at the same airmass 
as the target. The transformation equations are: 

H = h- (0.055 ± 0.005) x sec(z) + (23.430 ± 0.007) (C.l) 

with r.m.s.=0.040 mag for INGRID, 

H = h- (0.032 ± 0.003) X sec(z) + (24.727 ± 0.006) (C.2) 

with r.m.s.=0.01 1 mag for the Oct 2002 SOFI run, and 

H = h + (24.775 ± 0.009) (C.3) 

for the Aug 2003 SOFI run. Here the upper case indicates mag- 
nitudes in the standard system, and the lower case indicates the 
instrumental system. 



